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ABSTRACT 

We study the evolution of streams in a time-dependent spherical gravitational potential. Our goal is to establish what are the imprints 
of this time evolution on the properties of streams as well as their observability. To this end, we have performed a suite of test- 
particle experiments for a host system that doubles its mass during the integration time and for a variety of initial conditions. In 
these experiments we found that the most striking imprint is a misalignment of ~ 10° in the angular location of the apocentres of the 
streams compared to the static case (and to the orbit of the centre of mass), which only becomes apparent for sufficiently long streams. 
We have also developed an analytic model using action-angle variables which allows us to explain this behaviour and to identify the 
most important signature of time evolution, namely a difference in the slope defined by the distribution of particles along a stream in 
frequency and in angle space. Although a difference in slope can arise when the present-day potential is not correctly modelled, this 
shortcoming can be by-passed because in this case, streams are no longer straight lines in angle space, but depict a wiggly appearance 
and an implausible energy gradient. The difference in slope due to time evolution is small, typically ~ 10“^ and its amplitude depends 
on the growth rate of the potential, but nonetheless we find that it could be observable if accurate full-space information for nearby 
long streams is available. On the other hand, disregarding this effect may bias the determination of the present-day characteristics of 
the potential. 
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evolution 


1. Introduction 

Throughout the history of the Milky Way and in the context of 
the A cold dark matter cosmogony, many dwarf galaxies must 
have been d isrupted, leaving behind stellar streams ( [Helmi & 
White|[T999] ). Especially the Galactic halo probably contains 


galaxies ( jMartmez-Delgado et al.|2010[ Martin et al.|2014| ), most 
notably around M31 ( jlbata et al.|2001a l. 

On the modelling side, much effort has been put into under¬ 
standing the dynamics of streams and using these to infer the 
Galactic potential. Stream stars follow trajectories that are very 
similar to those of their progenitors ( jJin & Lynden-Bell||2007 


|Binney|[2008| ), albeit slightly offset. The evol ution of streams 
is relatively simple in action-angle coordinates ( [Helmi & Whitej 
[1999[[TTemaine|1999| ). Eurthermore, streams appear as distinct 


clumps in integrals of motion space ( [Helmi et al.[1999t [Helmi & 


many such relics of this assembly history ( [HelmipOOS] ). 

Tidal streams were predicted for the first time in the seminal 
work of [Toomre & Toomr^ ([1972[) who simulated the interac¬ 
tions of two galaxies. [Lynden-Beh & Lynden-Beh[ ( [1995[ ) and 
Johnston et al.[ ( [1996j ) put forward the idea of tidal streams lurk¬ 
ing in the Galactic halo, merely a few years after the discovery 
of the Sagittarius dSph by [Ibata et al.\ ( 199^ . Streams and de- 
bris from the Sagit t arius stream were found several years lat er by 
Ivezic et aL] ( [2000[ ); |Yanny et J^ ( [2000[ ); [Ibata et al.[ ( [2001b[), and 
a full-sky view of the Sagittarius stream was provided by [Majew^ 
ski et al.[([2003[). Other examples are the glob ular cluster streams 
([Grihmair et al.[[1995l, such as Palomar 5 ([Odenkirchen et al.[ 


[de Zeeuw[2000[ ), and they show the highest coherence in the true 
(underlying) potential. This can be used to constrain its char- 
acteri stic parameters ( [Penarrubia et al.[|2012[ [Sanderson et al. 
[2014[ ). The alignment of the angles and orbital frequencies can 
aid in determining the potential, as proposed by [Sanders & Bin- 
ney[ ( [2013a| . 


2001[ ) and NGC5466 ( Grihmair & Johnson|2006[ ). A wealth of 

new stre ams have been found in the Sloan Digital Sky Survey 
(see e.g. [Grillmair & Dionatos[[2006[ ), in the Tield of Streams’ 
( [Belokurov et al.[2006[). Streams have a lso been found with other 


The current cosmological model predicts significant evolu¬ 
tion in the mass content o f galaxies and of th eir dark matter 
haloes through cosmic time ( [Springel et al.[2005) ). This evolution 
may be directly measurable using stellar streams, given their sen- 
sitivity to the gravitational potential in which they are embedded 


([Johnston et al.| 1999} [Eyre & Binney|2009[ [Price- Whelan et al. 


2014t[Bonaca et al.[2014[ ). A first study was presented by [Penar¬ 
rubia et al.[ ( [2006[ ), who concluded that there are no discernible 


effects of evolution on the distribution of streams in the space of 
angular momenta and energy because they only reflect the poten¬ 
tial at the present day. However, [Gomez & Helmi[ ( [2010[ ) found 
that the structure of streams in frequency space does depict long- 
lasting signatures of time evolution. [Bonaca et al.[ ( [2014[ ) also 
found several biases when they attempted to derive the charac¬ 
teristic parameters of a time dependent potential using streams. 
They attributed these biases partly to the time evolution of the 
potential in their simulations. 

In the coming decade, the Gaia satellite ( [Perryman et al. 
[2001[ ), successfully launched in 2013, will provide an unprece¬ 
dented vast and detailed view of our Galaxy. With the correct 
understanding of streams, Gaia will allow us to address the as¬ 
sembly history of our Milky Way. This is the main motivation 
of this paper. Our goal is to establish what kinds of observable 
imprints remain on the stream properties and how to use them 
in recovering the evolutionary path of the Galactic potential. To 
this end we use the cosmologicahy motivated growth model of a 
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(2014) and follow the evolution 


of streams both numerically and using the action-angle formal¬ 
ism. 

This paper is organised as follows. In Sect. we explain 
the simulations of streams in a time-dependent potential. We 
analyse the simulations in Sect. and develop an analytic model 
based on action-angles in the adiabatic regime. This allows us 
to identify the signature of time evolution on a stream in Sect.|^ 
In Sect. we explore the prospects of observing this effect. We 
conclude in Sect.[6l 


2. Methodology 


To understand the behaviour of streams in time-dependent po¬ 
tentials, we take two complementary approaches. We first nu¬ 
merically study the evolution of groups of test particles that ini¬ 
tially resemble a satellite galaxy. We then attempt to model this 
evolution using a formalism based on action-angle variables. In 
Sect. |2.1| w e des cribe the set-up of our numerical experiments, 
and in Sect. |2.2| we briefly introduce the action-angle variables 
and their properties. 


2.1. Simulation set-up 

2.1.1. Evolution of the gravitational potential 



t(Gyr) 


Fig. 1: Increase of the enclosed mass for several orbits for our time- 
dependent potential with a growth factor of ag = 0.8. The dashed lines 
show the enclosed mass within the final value of their apocentre. Since 
the apocentric distances shrink with time, we show the mass enclosed 
within the time-dependent apocentre values with solid lines. The inte¬ 
gration time span for the streams is ^ ~ 7 - 8 Gyr, and the enclosed mass 
approximately doubles over this timescale. 


We modelled the evolution of a time -dependent Navarro- 
Frenk-White potential ([Navarro et al.|1997~l N FW) using the pre¬ 
scription models from jBuist & Helmi p014| ). This model gives 
rise to an inside-out smooth growth that well reproduces the 
average mass accretion history of the Milky-Way-sized haloes 
from the Aquarius simulations (jSpringel et al.|2008l|Wang et al.j 
[MTT] ). An important characteristic of this model is that mass 
growth is positive at each radius (i.e. there is no inward flow of 
mass between neighbouring mass shells). The scale radius and 
the scale mass Ms (mass enclosed within r^) of such a halo vary 
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according to 


Ms(z) = Ms,o exp (-2agz), 


rs(t) = rs,o 


Ms(t) 


s,0 


1/r 


(1) 


where z is the redshift, and Ms,o and rs,o the scale mass and scale 
radius at the final time (z = 0, in a cosmology with Qm = 0-29 
and /zioo = 0.71). We used y = 2 to grow inside-out for an NFW 
potential (jBuist & Helmi||2014|). The parameter Ug determines 
the formation epoch, and higher values correspond to more re¬ 
cent formation. We here explore a range of values of Ug up to 
Gg = 0.8, which would represent a quite recent formation epoch. 

We modelled the evolution of streams in this time-dependent 
potential, but we also studied their behaviour in the static case. 
To make a fair comparison, we fixed the scale mass and scale 
radius at the final time to Ms = 5 x 10^^ Mq and = 12 kpc for 
all our experiments, independent of their time evolution. 


2.1.2. Orbits 


Table 1: Orbital properties at the final time 


Orbit 

^apo 

(kpc) 

^peri 

(kpc) 

T/Tjnax 

Qr 

(rad/Gyr) 

(rad/Gyr) 

Q(f,/Qr 

01 

90 

5 

0.22 

11.38 

6.96 

0.61 

02 

70 

5 

0.26 

15.06 

9.25 

0.61 

03 

60 

5 

0.29 

17.79 

10.96 

0.62 

04 

50 

5 

0.33 

21.54 

13.32 

0.62 

05 

90 

10 

0.39 

11.00 

7.25 

0.66 

06 

40 

10 

0.64 

24.97 

16.60 

0.66 

07 

70 

15 

0.61 

13.69 

9.45 

0.69 

08 

30 

20 

0.96 

25.89 

18.15 

0.70 

09 

90 

25 

0.71 

9.68 

7.01 

0.72 

OlO 

40 

30 

0.98 

17.83 

12.93 

0.72 

on 

60 

30 

0.90 

13.21 

9.68 

0.73 

012 

90 

30 

0.78 

9.25 

6.81 

0.74 

013 

50 

40 

0.99 

13.32 

9.89 

0.74 

014 

80 

40 

0.90 

9.36 

7.03 

0.75 

015 

80 

50 

0.95 

8.53 

6.49 

0.76 


Table 2: Progenitor properties 



a-pos (pc) 

o-vei (kpc/Gyr) 

‘Carina ’ 

100 

5 

‘Sculptor' 

300 

10 

‘Sagittarius' 

700 

25 


We placed the satellites on 15 different orbits as listed in Table 
[2 We chose this number to have a variety in distances, orbital 
timescales (and hence the degree of adiabaticity), and eccentric¬ 
ity of the streams. The apo- and pericentres at the final time were 
chosen such that the particles spend most of their time in the 
outer halo and beyond the scale radius of the potential. For a fair 
comparison of each experiment the final 6D position of the satel¬ 
lite’s centre of mass is the same for the static and time-dependent 
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cases. We then integrated the orbits backwards in time for about 
7-8 Gyr and the more circular orbits for 6-7 Gyr. This range was 
chosen to ensure that all orbits are bound at early times in the 
time-dependent case. 

We chose to start our orbital integrations when the satellite’s 
centre of mass was at its first pericentre to represent the first in¬ 
teraction with the Galaxy. This implies that the various orbits 
have slightly different total integration times, but these only dif¬ 
fer up to a maximum of 0.5 Gyr. Within this integration time, the 
enclosed mass within a fixed radius increases on average with 
more than a factor 2 in our halo, as seen in Fig. The actual 
orbits experience a somewhat smaller increase in enclosed mass 
because they respond to the halo growth by shrinking slowly. 


2.1.3. Progenitors 

We distributed the particles in the satellites assuming they follow 
an isotropic Gaussian in position and in velocity space, charac¬ 
terised by dispersions CTpos and CTyei. We used 10,000 particles 
centred on the centre of mass of the satellite, which was placed 
on the orbits described previously. We did not include self¬ 
gravity in our simulations. We considered three different progen¬ 
itors, which we called ‘Carina’, ‘Sculptor’ and ‘Sagittarius’ be¬ 
cause they have properties reminiscent of these dwarf spheroidal 
satellites of the Milky Way (see e.g. [Martin et al.||2008l |Wolf| 
jet al.|2010| see Tablel^. 


2.2. Action-angle coordinates 
2.2.1. Generalities 


A particularly useful description of the evoluti on of streams may 
be obtained using action-angle coordinates ( |Goldstein||1950 


[Helmi & White] 1999} [Tremaine] 1999[ ). In action-angle coordi¬ 
nates, the actions are the momenta and the angles are the coordi¬ 
nates. The actions are integrals of motion that are adiabatic, or in 
other words, invariant under slow changes of the host potential. 

In this section, we only work with a time-independent poten¬ 
tial. For a single star, the actions can be found from 



( 2 ) 


with Pi the conjugate momentum of coordinate qt, and ^ to in¬ 
dicate that we integrate over a full period of coordinate qi. In a 
spherically symmetric system the actions are the radial, latitudi¬ 
nal and azimuthal actions (/^, Jq, and respectively), and the 
angles are 6r, Og, and 6(p, and represent the phase of the orbit in 
the r, 6, and 0 coordinates, respectively. 

For only a few potentials algebraic expressions of the actions 
can be derived, such as for the isochrone potential (of which 
the Kepler potential and harmonic oscillator are limiting cases). 
Typically numerical methods are therefore required to find the 
actions ( [Binney & Tremaine[[2008| ). Furthermore, actions are 
only directly available when the Hamilton-Jacobi equations are 
separable, as in the case of spherical and Staeckel potentials 
( [Goldstein|1950[[de Zeeuw|1985[ ). In more general axisymmet- 
ric and triaxial potentials, approximations to the actions and an- 
gles must be used (see e.g.[Kaasalainen & Binney 1 19941 Binney[ 


& McMillan|20Tll[Binney|2012|[Sanders & Binney|2014[[2015t 

Bovy|2014]K 



0 2.7T 47T 

e^(\) (rad) 

Fig. 2: Relation between the coordinates and the angles in the plane of 
the orbit. The top panel shows that the radial angle 6r increases by In in 
one radial oscillation. In the bottom panel we have plotted 0 (red) and 
6(p (blue) as functions of 6r. After one radial period, both angles reach 
the same value 2nQ(f,IQr- 


Hamilton’s equations for the // and 6i are 


Ji 

Oi 


dH(J) 

dOi 


= 0 , 


dHQ) 

dJi 




( 3 ) 


and H = H{i) only because the actions are invariant. The Q/ 
are constants that correspond to the orbital frequencies and also 
only depend on the //. The resulting equations of motion for the 
Oi are 


6i — ^it ^/(O), 


(4) 


with OiiOi) the initial phase of i-th coordinate. 

To determine the angles, we used a canonical transformation 
of the second kind W(q, J) ( [Goldstein[l 95 Oj [Binney & Tremaine 
[2008[ ), such that 


Pi 


dWjqJ) 

dqi 


dW(qJ) 

dJi 


(5) 


The generating function for the transformation in a spherical po¬ 
tential is given by (adapted from [Binney & Tremaine|200^ 


W(q, J) = W^(0, J) + Wdd, J) + W,(r, J) 

d(l)Pcp(Jcp) + I d§ pdJcp, J&) + I dr Pr(Jr, Jcp, 


( 6 ) 


where the integration is along the trajectory of a particle in phase 
space. We note that the three parts of the generating function are 
indefinite versions of the action integrals without a factor of 271, 
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Fig. 3: Distribution in action (left), frequency (second panel) and angle space (right) for particles in a stream, and centred on the centre of mass 
of the progenitor system. The colours represent the energy gradient, where red is for the trailing and green for the leading arms. 


and each increases therefore by InJt in their corresponding peri¬ 
ods This equation has an algebraic expression in the case 

of an isochrone potential and has to be solved numerically oth¬ 
erwise. For the full integrals needed to compute 6i in a spherical 
potential we refer to Appendix [A| 

In Fig. we show an example of how the radial and az¬ 
imuthal angles (Or and O^p) are related to the radial and azimuthal 
coordinates (r and 0) in the orbital plane. Since = 0/^ and 
using Eqs. ^ and ( [bf] 


0(p 


dW dWr 

— = 0 -r-. 

dJ^ ^ dJ^ 


(7) 


The term has a non-secular oscillation with frequency 
oJs 


and vanishes at apo- and pericentre (Or = {0, ±n}). The secular 
behaviour of 0 and 0(p with time is therefore the same, as can also 
be seen in Fig.[^ 


2.2.2. Streams in action-angle coordinates: 
time-independent case 

A stream is created when particles drift away from the progen¬ 
itor, which means we should look at the relative phase with re¬ 
spect to the progenitor centre of mass. For example, the relative 
azimuthal phase of the k-th particle is 

^(0 - ^”(0 = A0j(r) = A0^(O) + 

«AQ^r, (8) 

where A0^(0) is the separation in angles at the initial time, and 
where the second line corresponds to the case in which the pro¬ 
genitor is small and the stream has evolved for a sufficiently long 
time. 

The AQ^ = depends on the distribution in ac¬ 

tion space such that a broad action spread results in a wide 
range of frequencies. This implies that larger progenitors gen¬ 
erate longer streams than smaller ones if integrated for the same 
amount of time. 

In Fig. we show the structure of a stream in action-angle 
coordinates for the time-independent case. The colours represent 
the energy gradient, such that green represents the leading tail 
and red the trailing tail. The left panel shows the action distribu¬ 
tion, which depends both on the initial conditions (Xpos and CTyei 

dWe 

^ —— = 0 on the orbital plane. 
dJ(p 
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and on the initial orbi{3 This distribution also remains invariant 
in time in an adiabatically evolving potential. The second-left 
panel of shows the structure in frequency space, while the third 
corresponds to angle space. Since angles are In periodic vari¬ 
ables, in the right most panel we show the angles modulo In. 
The different streaks signify the number of radial and azimuthal 
wraps the stream has. We note that the slope of the lines in fre¬ 
quency and in angle space are the same in this time-independent 
potential ( [Sanders & Binney]|2013a| ) when the initial spread in 
angles can be neglectec0 


AOfj) 

Afr 


AQ.(p 

AQr 


= cst. 


(9) 


This directly shows why all particles are distributed along 
straight lines with the same slope in angle and in frequency 
space. 


3. Analysis of the test-particle simulations 

The results of the simulations for the three progenitors on the 
various orbits for the time-dependent potential with growth fac¬ 
tor = 0.8 are shown in Fig. The colour coding indicates 
the leading (green) and trailing (red) arms of the stream. For 
comparison we have over-plotted the orbit of the progenitor in¬ 
tegrated forward and backwards in time in the present-day poten¬ 
tial. The time of integration for that orbit is between 1 to 3 Gyr 
and is chosen to roughly match the length of the corresponding 
stream. Note that this is much shorter than the time of integra¬ 
tion of the particles that form the stream, which is approximately 
8 Gyr. 

Figure clearly shows that larger progenitors give rise to 
longer and wider streams. The smaller (Carina) progenitor gen¬ 
erates thin, short streams, which makes the separate loops of 
the stream easily recognizable by eye, while the streams from 
the larger progenitors, depending on the specific orbit consid¬ 
ered, are wider and longer, much more phase mixed, and some¬ 
times difficult to discern. As we show below, a clear imprint 
of the time evolution of the potential is apparent in the angular 
location and in the variation of the maximum distance of each 
‘petal’ of a stream (the latter being related to both the shrink¬ 
ing of the progenitor’s orbit because of mass growth and to the 
energy gradient along a stream). This effect is more easily dis¬ 
cernible for heavier progenitors since they produce more loops, 

^ Although here we plot L = -h l/^l, but in the orbital plane = 0. 

^ This equation holds for every particle that is released at the same time 
from the progenitor. 
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Fig. 4: Streams from three progenitor sizes for each of the 15 orbits evolved in the time-dependent ag = 0.8 potential. The left, middle, and right 
panels correspond to the ‘Carina’, ‘Sculptor’, and ‘Sagittarius’-like progenitors, respectively (see Tablej^for their properties). Particles in green 
are more bound to the host galaxy (leading arm), while those in red are less bound (trailing arm) than the progenitor. The dashed curves are the 
progenitor orbits evolved in today’s potential, with the green integrated forward in time for the leading arm and red backwards for the trailing arm. 
The current position of the progenitor is indicated with a red cross. 
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Fig. 5: Radial vs angular distribution of particles in the orbital plane for experiment 07 for the different progenitors in the time-dependent case. 
The dashed curves show the orbit of the centre of mass of the system integrated long enough to roughly reproduce the lengths of the streams. The 
variations in the apocentric distances are a reflection both of time evolution and of the energy gradient present along a stream. The difference in 
angular location of the apocentre, on the other hand, is a clear imprint of time evolution, as can be seen by comparing to e.g. Fig.|^ 


that is, longer streams, as can be seen in Fig. On the other 
hand, thinner streams have the advantage that they depict their 
mean orbit much more clearly, which is the reason we mainly 
focus on the ‘Carina’ progenitor in what follows. 

The differences in the evolution for time-dependent and 
static potentials are shown in Figs. and [7] for the Carina-like 
progenitor, depicted as before on the orbital plane of the progen¬ 
itor. We note that in all cases the stream in the time-dependent 
potential is longer than in the static case, as can be seen by com¬ 
paring the left and middle panels. 

The length of a stream depends both on the orbital period 
and on the initial extent in phase- and configuration space of the 
progenitor. In the time-dependent potential the enclosed mass at 
initial times is lower, and the progenitor starts farther out. This 
results in lower orbital frequencies, or longer periods. Since the 
rate of divergence of particles in a s tream in a spherical pot ential 
is proportional to {tlP)^ oc ( Helmi & White|p^^ ), this 

would imply that streams should actually be shorter for a fixed 
integration time t. However, in a time-dependent potential the 
spread in actions (or integrals of motion such as the energy) is 
broader for given initial CTpos and (Tyei (see Fig.[^. 

For an initially ‘compact’ ensemble of particles, we may 
use a linear transformation between action-angle coordinates and 
Cartesian coordinates to derive the initial spread in action-angle 
space: 


^e 

AJ 


de 

de 

5q 


dl 


d<\ 



Aq 

Ap 


= T 


Aq 

Ap 


( 10 ) 


where the elements of the transformation matrix T will depend 
on second derivatives of the generating function evaluated at the 
centre of mass of the progenitor ( [Helmi & White] 1999] ), as de¬ 
scribed in detail in Appendix We see that many of the terms 
depend on which is larger in a shallower potential and 

farther out in the potential, confirming that the initial spread in 
actions is greater. Therefore the particles in the progenitor inte¬ 
grated in the evolving potential can spread more and give rise to 
longer and wider streams. 
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In our initial set-up we took the progenitors to have the same 
CTpos and CTyei. However, their initial distribution in action-angle 
space is different for the static and time-dependent potentials, 
as discussed above. This implies that their final distribution in 
action (or energy) space will also be different. This is shown for 
the progenitor placed on two of our orbits in Fig. Here we 
have plotted the distribution at the final time of AA and AL in 
the top panels for the static (left) and time-dependent (middle) 
potentials for the ‘Carina’-like progenitor. 

An alternative is to consider that the particles have the same 
initial AA and AL, whether they are evolved in a static or time- 
dependent gravitational potential. To establish the effect of such 
a change in the initial configuration, we set up progenitors with 
these properties and evolved them in the static potential. The 
right panels of Fig. show the distribution of particles evolved 
in the static potential, but set up with the same initial conditions 
in action-angle space as in the time-dependent case. If all orbits 
were in the adiabatic regime, the middle and right panels should 
give identical results. However, some of the orbits are not, which 
can be seen in the streams from orbits Oil and 015, where the 
action distributions are not exactly the same. 

The streams resulting from these different initial conditions 
evolved in the time-independent potential are shown in the right 
panels of Fig.j^for each of the orbits. As anticipated, the streams 
in the time-dependent potential (middle panels) are now shorter 
than the streams in the time-independent potential starting from 
the new initial conditions. This characteristic is now solely due 
to the evolution of the potential. 

An interesting feature visible in Figs. andthat was men¬ 
tioned earlier is that the stream-orbit misalignment differs in 
the static and evolving cases. For example, as can be seen for 
experiment 07, the mean orbit traces the stream in the time- 
independent case relatively well. On the other hand, it is clear 
that the angular position of the rosette petal of the stream and the 
progenitor orbit at one radial period behind (or ahead) are offset 
from each other. This is apparent in all experiments, to a lesser 
or greater degree. This offset is systematic, can be as large as 
10 degrees, and is a first direct indication of the effect of time 
evolution on the gravitational potential. 
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Fig. 6: Streams from the Carina progenitor in the evolving and in the non-evolving potential. The streams in the left and centre have the same 
initial distribution in configuration and velocity space. Those in the right panels have the same initial action-angle distribution as the time- 
evolving potential (middle panel), but are evolved in a static potential. The colour coding is the same as in the previous figure. The stream-orbit 
misalignment can be seen by looking at the azimuthal angle of the petals starting from one radial period away from the progenitor, whose position 
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Fig. 7: Radial vs angular distribution of particles in the orbital plane for experiment 07 for the Carina-like progenitor in the static case (left), 
time-dependent case (middle), and static case with the same initial action-angle distribution as in the middle panel. The dashed curves show the 
orbit of the centre of mass of the system integrated in the present-day potential for the left and right panels, and in the evolving potential for the 
middle panel. For the static cases, the variations in the apocentric distances are a reflection of the energy gradient present along the stream, while 
a second effect is present in the middle panel because of time evolution. In this case, the difference in angular location of the apocentre is also 
much more pronounced and systematically increases for older wraps. 
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Fig. 8: Action-space distributions for a sample of four streams from Fig.|^ measured with respect to the centre of mass of the progenitor at the 
final time. The colours of the particles represent the energy gradient, such that green is the leading arm and red the trailing arm. In the static cases, 
the initial and final distributions coincide exactly, but in the time-dependent case there can be small deviations if the orbit is not in the adiabatic 
regime, as can be seen by comparing the middle and right panels for orbits Oil and 015. The left panels show a smaller spread because the 
time-dependent cases start farther outwards and in an initially much lighter potential. 


4. Analytic models 

4.1. Action-angle coordinates in an adiabatically evolving 
potential 

In this section our aim is to extend the use of action-angles to 
an adiabatically changing spherical potential. The action-angles 
are still a valid canonical coordinate system in a time-dependent 
system, but the equations of motion in action-angle coordinates 
are more complicated, as we show below. 

In the models we have considered so far, the gravitational 
potential is made time dependent by making its characteristic 
parameters, such as mass and scale, a function of time. Its overall 
shape or functional form remain the same. Therefore angular 


momentum is still conserved because of the assumed spherical 
symmetry, but the radial action Jr may vary with time. 

For a time-dependent potential, the generating function that 
allows the transformation between Cartesian and action-angle 
spaces is an explicit function of time 

W(qJj) = W(qJ,a(t)), (11) 

where the parameters of the potential are in the vector a. The ap¬ 
propriate Hamiltonian in action-angle coordinates then becomes 


, dW dW 

= H(J,a)+ — (qJ,a) = H(J,a)+a—(qJ,aX (12) 
dt da 

which is the original Hamiltonian perturbed by the partial time 
derivative of the generating function ( |Goldstein||1950[ section 
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11-7). This extra term naturally vanishes in the time-independent 
case. The new equations of motion are 


SH' _ a atr 

. dH' __ ^ dW, ^ ^ 

^i = -^ = or) + or— —(q, J, or). 
oJi oJi oa 


(13) 


Compared to the original equations of motion given in Eq. 
there are a few correction terms that depend on the rate of change 
of the characteristic parameters of the potential a. 

Since in general we are interested in the mean increase in the 
actions and angles, we take an average over one period, where 
we assume that a changes little over one period ( |Goldstein|195Q 
I VandervoortI 1961 1 ) 

/.X \ r d dW^ ^ 

~~¥ f -^^(^^3,a)dt + 0{a^,a). (14) 

Ti Jt, d0i da 


The three parts of the generating function {Wr, Wg and W^) in¬ 
crease by 2nJi in the corresponding period Ti. This means that 
each of the Wi is a periodic function, and we can do a Fourier 
expansion 


k 

However, the whole first term in Eq. ( p^ vanishes since 

f Zj + 0(dj, aj 

= 0(d^, a). 


(15) 


(16) 


where we use the Einstein summation convention for the param¬ 
eters am- The actions are therefore invariant up to second order. 
|Vandervoort| \\96\) showed that the actual condition for adia¬ 
batic evolution for each of the am and each of the orbital periods 
Ti is 


— Ti « 1, (17) 

which states that the timescale of change should be much longer 
than the orbital period. 

On short timescales, Jr shows a periodic oscillation that de¬ 
creases in amplitude when the changes in the parameters become 
smaller, as shown in Fig. for one of our orbits. The final time 
value of Jr returns to the initial value, indicating we are in the 
adiabatic regim^ In the case of non-adiabaticity, full integra¬ 
tion of Eq. •ini' is needed to model the behaviour of the system 
in action-angle space. 

In an analogous way as for the actions, we compute the effect 
for the mean change in the angles 

{e) ^ (Q,(J, «)> + J’ «)’ ( 18 ) 


where we can exchange the order of the derivatives 


d dW d dW dOi 


(19) 


We note that if the particles were not released from pericentre, we 
would notice the initial variations in the action. 



Fig. 9: Response of the radial action and the radial frequency to a 
time-dependent potential for orbit 07 (Table [^. For this orbit, the ac¬ 
tion shows a non-secular oscillation that decreases in amplitude as the 
rate of change of the potential decreases. The value of changes sig¬ 
nificantly with time, and here there are also oscillations that become 
smaller as the change rate of the potential decreases. The period of the 
oscillations corresponds with the radial period. The evolution of A has 
been quite adiabatic because it oscillates around a limiting value, which 
is the starting value when the particles were released at pericentre. 


The result is that only the term dependent on the (evolving) fre¬ 
quency remains 

{4)~(Q,(J,a)>+<9(a^a). (20) 

The angles as a function of time may be computed from 

^,(0^^/(0)+ r ^i{i,a{t))du (21) 

Jo 

where like in the time-independent case, ^/(O) is the phase at 
the initial time. The frequencies Q/(0 are understood to be the 
instantaneous frequencies at every time step. Equivalently, we 
can numerically integrate the orbit in Cartesian coordinates and 
find 6i(t) from the instantaneous coordinate transformation at 
all times t using the generating function. The advantage of the 
action-angle description is that it stresses that the phase angles 
(modulo 2;r) depend on how the frequencies have changed with 
time. 

4.2. Streams in action-angle space 

4.2.1. Sensitivity to the gravitational potential 

We have mentioned that if the initial spread in angles is negligi¬ 
ble, then the slopes in the angle and frequency spaces should be 
equal for a static potential ( [Sanders & Binney|2013b|a| ), as can 
be seen from Eq. 0, if computed in the true potential. An ex¬ 
ample of the distortion in angle space for a stream now evolved 
in a time-dependent potential where the angles (and the energy) 
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Fig. 10: Angle space of stream 07 (see Table[^ evolved in a time-dependent potential, but where the angles have been computed with different 
(incorrect) parameters. The central panel shows the correct parameters, while is changed vertically and Ms horizontally by 10% and 30%. The 
colour coding indicates the energies, with green representing the leading arm (more bound) and red the trailing arm (less bound). 


were computed using incorrect parameters is shown in Fig. 10 


for experiment 07. The central panel corresponds to the pa¬ 
rameters in the true (final) potential, while we varied in the 
vertical and Ms in the horizontal direction. The parameters in 
the eight panels around the central panel were changed by only 
10%, and even in this case, some deviations from the expected 
straight lines are visible. When the parameters were changed by 
as much as 30%, as in the outer panels, the deviations become 
much more pronounced and are very strong. We note that similar 
results are found when the stream is evolved in a static potential. 


The more incorrect the potential, the more we see wiggles 
in angle space. Similar distortions are also present in the fre¬ 
quency spaces and action spaces. Furthermore, the energy gra¬ 
dient along the stream is no longer preserved as the parameters 
of the potential are varied, which implies that the most bound 
(least bound) particles in the specific trial potential are no longer 
‘found’ at the end of the leading tail (trailing tail). Although the 
streams have a normal appearance in physical space, the coordi¬ 
nate transformation to action-angle space and the computation of 
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the energy are incorrect when the wrong gravitational potential 
is assumed. This leads to the distorted appearance and broken 
energy gradient seen in Fig.[T^ Clearly, these imprints will only 
be obvious for sufficiently long streams, that is, those more ex¬ 
tended than 271 in angle space. 


We note that the distortions in angle space are less pro¬ 
nounced for the panels located along the diagonal that runs from 
bottom left to top right in Fig. 10 for which the enclosed mass 
within the orbit is the same for the given and r^. 


The degree to which these distortions and broken energy gra¬ 
dients manifest also depends on the type of orbit the stream pro¬ 
genitor has followed. This is shown in Fig.[TT]for a scale radius 
that is 70% of the true value and for different scale masses M^. 
The top panels correspond to experiment Oil, which is on a rel¬ 
atively circular orbit of apocentre-to-pericentre ratio of 2, while 
the experiment in the bottom 01 has a relatively radial orbit, with 
apocentre-to-pericentre ratio of 18. 
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Fig. 11: Angle space of streams Oil (top) and stream 01 (bottom) evolved in a time-dependent potential, but where the angles have been 
computed with different (incorrect) parameters. We note that although the distortions in angle space are less pronounced for the more radial orbit 
01, the energy gradient (indicated by the colour gradient and coding used before) is not preserved. 


4.2.2. Sensitivity to time-(depen(dence 


In the time-dependent potential the exact correspondence be¬ 
tween the slopes of the spreads in the frequencies and in the 
angles is broken, meaning that Eq. (|^ is no longer valid. This is 
because the evolution of A^^ and A6r depends on integrals over 
time, that is, 


A6(plA6r ~ 


J" A^^dt I 


AQ^rdt 


( 22 ) 


(assuming the initial angle spreads are negligible). 

In Fig. we compare the lines fitted to the frequency and 
angle spaces for the stream from orbit 07 in the time-dependent 
and static cases for the correct final potential. The top panel 
shows that for the time-dependent case, the initial slope in fre¬ 
quency space (grey line) is steeper than at the end of the simu¬ 
lation (blue line). This is because the orbit shrinks as the mass 
increase^ In the time-independent case we see, as expected, 
that both lines coincide. Another difference between the static 
and evolving case is that the stream-orbit misalignment changes 
(compare the green and red dashed lines), which is only caused 
by the changes in the angle space due to time evolution. The 
behaviour in angle space (indicated by the red dashed line) may 
be understood from the fact that the angles’ slope is like a time 
average of the frequencies’ slope, and is therefore expected to 
lie between the initial and the final frequency slope as observed. 
Therefore, a clear signature of time evolution is a difference in 
the slope derived for the angles and that derived for the frequen¬ 
cies, even if they are computed using the present-day gravita¬ 
tional potential. 

^ We recall that the ratio of angular to radial frequencies for any gravi¬ 
tational potential is limited by the homogeneous sphere case, for which 
= 1/2 and the Kepler case, for which = 1- As mass 

increases and the orbit shrinks, we may say, effectively, that the stream 
moves farther away from experiencing a Kepler potential (and hence a 
steeper slope) and closer to the homoge neous sphere (a shallower slope, 
as also seen in [Gomez & Helmi|2010| ), although the NFW potential is 
neither of these limiting cases. 


(rad/Gyr) 

0 2 4 6 



Fig. 12: Lines fitted to the angle space (red dashed line) and the fre¬ 
quency space (blue line) of stream 07 evolved in a time-dependent (top) 
and static potential (bottom) using the same (correct) final potential. 
Additionally, the green line shows the track of the progenitor assuming 
the potential is static (i.e. 0cm(O = and the grey line shows 

the line fitted to the frequency space at the beginning of the simulation. 
The inset labels show the numerical values of the aperture angle (mis¬ 
alignment) between the stream and the (static) progenitor orbit 'Fs,o. 
and between the stream in angle and in frequency space (see 

Eqs. ( [^ and ([^). 
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Fig. 13: Hessians Hij, time-averaged Hessians Qij (both in units of 10"^ kpc"^), and their ratios hi and qi for orbit 07. Left panel: the dash-dotted 
lines show Hij in the static potential, the solid lines the Hij in the time-dependent potential, and the dashed lines Qtj in the time-dependent potential. 
We note that |7/^r| > and similarly for Qij. The middle panel shows that at all times and that both increase with time. On the 

other hand, hr < qr at all t, and both decrease with time. Right panel: the ranking of the hi for the different orbits. Both h^p and hr are ranked by 
the ratio 0^/0^, indicating the dependence on the type of orbit. The qi are ranked in the same way. 


We now explore an analytic model to describe the behaviour 
of the frequency and angle slopes more generally. We may de¬ 
rive the slope in frequency space S (AQ) by making a Taylor ex¬ 
pansion in the actions near the centre of mass of the progenitor 
( Helmi & White] 1999| ) 


AQf = ^Ay) + 0(Ayf). 


(23) 


The expansion is done for the ^-th particle, all terms are evalu¬ 
ated at time t. The derivative is evaluated for the progenitor orbit 
and we used the Einstein summation convention on the subscript 
indices. We recall that 


d^i _ d^H _ ^ 

Wj ~ dJidJj ^ 

Working out the frequency spreads, we find 


(24) 


AQf(0 « Hrril)A4 + HremA + 

AQ^(0 « H^rmJr + + H^ 41 )AJ'^. 

The cross terms are second derivatives of the Hamiltonian and 
therefore symmetric, Hij = Hji. Furthermore, in a spherical po¬ 
tential, the Hamiltonian is only a function of E and L, so that for 
The evolution of the slope of S (AQ) can then 

be found from 


A(AQ) = 


AQ^(f) 

AQ*(0 


where /?j = 
h^{t) = 


A4 

AJl + AJl ^ 

^ ^ HrXt) Aj] 

A/^ + AJ’I 

a4 

A^ + Ay^’ 

H^rity 


hr{t) = 


Hrrit) 


+ h^{t) 

l+R'lhrity 


(25) 


^ This is only true if > 0, otherwise every derivative w.r.t. incor- 
porates a term signi/^) because L-jQ + \Jj)\. Throughout this work we 
assume 7^ > 0. 
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Here /?j is the ratio of the actions and hi the ratio of the Hes¬ 
sians. The Hessians Hij and the Hessian ratios hi are shown 
in Fig. and are ranked as But most 

importaWy, the resulting ratio hr increases with time, while 
decreases with time. Assuming that the action ratio remains 
constant (adiabatic limit), we find from this equation that the 
slope in frequency space S (AQ) decreases with time. 

For the angles of the ^-th particle we can use the following 
approximation 

A6»f(0 = A6»f(0) + r An'l(t)dt 
Jo 

= A0f(O) + Hij{t)Aj){f)dt + 0{Ajf) 
t)Aj]{t)dt 

Hij{t)dt, (26) 

where we have neglected the initial angle spread and in the last 
approximation assumed adiabaticity of the actions. Therefore 
A/^ may be evaluated at any time step, although we typically 
use the values at the final time in our computations. 

We now focus on the slope in angle space. If we define 
y Hijdt = Qij (the time-averaged Hessian), then 



A0^(r) « {Qrr(t)A4 + QrB{t)A4 + Qr^Aj'^^ U 
A©^(r) ^ {Q^r{t)A 4 + Q^mfe + Qpp{t)Af^t. 
The slope in the angles can then be found from 


s{^e) = 


A00(O 

A0,(O 


where ^^(0 = 


+ ^^(0 
i+?r(o/?r 
Q<p<p{t) 

Qprit) ’ 


qr{l) 


Qrrit) 

Qprit)' 


(27) 


Here the qi are the ratios of the Qij. The ranking of the Qij and 
the qi behaves very similarly to that of the Hij and hi, as shown 
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in Fig. In the same manner as for the frequency slope, we 
infer that S (A6) is also a decreasing function of time. 

We are now ready to compute the difference between S (A6) 
and S (AQ) by using that the difference between the qi and the hi 
is small: 


qr{t) = hr{t) - Erit), (28) 

^ 0(0 = h(j){t) + 6^(0, (29) 

where the appropriate sign was chosen such that the 6/ are always 
positive. We define for computational ease 

We may expand the inverse of the angle slope T{M) 


1 + R]qr 1 + R]hr - R\er 
R\ + ^0 R\ + hip + 6^ 

1 + R'^jhr - R]er ( R'l + h^ \ 

Rj+hip \Rj + hp + ep) 

l+R’lhr-R'jer ( _ £0 ^ 

R'i + h^ ( R) + h^] 

1 + R]hr - R\er 1 + R'lhr 
R!l + h^ ^\R] + h^f' 


(31) 


Using Eq. ( [25) ) to compute T(A6), we can derive the difference 
of the inverse slopes 


T(A6) - r(AQ) 


S(AD.)-S(A6) 

S(A6)S(An) 

R]er 1 + R]hr 

^ + hcj) 

T{A^l)e^+R]er 

-^- — < 0 . 

R^ + hcj) 


(32) 


All the terms in the fraction are positive, and r(AQ) is there¬ 
fore always greater than T(A6). This means that S (A6) is always 
larger than S (AQ), and this difference will increase with time. 
The longer the system is evolved in a time-dependent potential, 
the larger the difference in slopes becomes with respect to the 
same final static potential. This also means that streams that fell 
in at early times will show a stronger signature of the separation 
between the angle and frequency relations than streams that fell 
in more recently. 

In Fig. we show the predictions of the slope differences 
for the angle and frequency spaces for the orbits in TableWith 
increasing amount of time evolution by changing the growth pa¬ 
rameter Qg, the difference in slope also becomes larger. Another 
interesting prediction is that the more circular orbits have a larger 
slope difference, although this is at a small level. We gauge 
the magnitude of the angle-frequency misalignment by using the 
aperture angle ao between the fitted lines in the angle and 
frequency space (as also used by [Sanders & Binn^ ( |2013b]a| )) 


AO = 'Fa^ - 'Fao ^ tan ('Fa^ - 'Fao) = 


S(A6)-S(A^) 

l-^S(A6)S(Any 

(33) 



Fig. 14: Slope differences between the streams in angle and in fre¬ 
quency space as a function of growth parameter as predicted by our 
analytic model using Eq. and p6^ for our experiments. 


where S (i) = tan('F/), and we have used that the difference in 
slopes is smalQ With this definition of the angle-frequency mis¬ 
alignment, typical values of'F a6»,ao are around 1.5° for ag = 0.8. 

It has been discussed earlier that a stream does not exactly 
follow an orbit, and the small difference between these two tra¬ 
jectories is known as the stream-orbit misalignment. The char¬ 
acteristic magnitude of this can be found from 

= = ( 34 ) 


where S (^cm) = tan 'F^^^ is the slope of the straight line traced by 
the centre of the mass of the progenitor in angle space for a static 
potential (see also Fig. [T^. The typical order of magnitude of 
'Fs, o for the Carina progenitor is about 1 degree, similar to what 
was found in [Sanders & Binne'y| ( |2013b| ) for a static logarithmic 
axisymmetric potential, and also similar to values of the stream- 
orbit misalignment in the isochrone potential of [Eyre & Binney 
( [201 1| ). We note that the angle-frequency differences in slope 
that we find are of the same order of magnitude as the stream- 
orbit misalignment. 

When comparing simulations with and without time evolu¬ 
tion, the stream-orbit misalignment can also be used as an in¬ 
dicator of time evolution. In Fig. we saw that the angular 
position of the stream and progenitor orbit petals one radial pe¬ 
riod away from the current progenitor location are different for 
the time-dependent case. This angular separation can be de¬ 
rived as follows. For the progenitor, the location of the petal 
is simply the precession of the orbit in one radial period, that is, 
i/^cm = = 27TQ(j)IQr’ Fof thc Stream, this difference in 

radial angle from the progenitor is A0^ = In, and therefore the 
angle of the petal of the stream is j/^stream = 27iS(A6). Therefore 
the stream-orbit misalignment in azimuthal angle is 


A^S.O ~ lAstream “ lAcm = 27r (S (A6») - S (6»cm)) ■ (35) 

Increasingly higher growth rates of the potential will yield also 
larger angular stream-orbit separations because S (A6) changes 

^ In these equations we have ignored the slope in the r-i^ angle and 
frequency spaces, because the slopes are very close to those of the r-0 
spaces. 
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Fig. 15: Examples of the slope differences between the angle and frequency distributions (top panels) and histograms of the spread of 
A^^(O)/A0^(t/) (bottom panels). In the top panels, the blue crosses indicate values from the full simulation, while the green lines are the pre¬ 
dictions from our analytic model. The red crosses are the slopes measured when the initial spreads from the simulations are removed. All orbits 
now align with the model when ag = 0. For large ag the more circular orbits such as orbit 015 deviate strongly from the model because of 
non-adiabaticity. 


by the amount in Fig.[^(5' (AQ) and S (0cm) are by construction 
the same at the final time). By multiplying this by 360° to find 
the difference at the first petals next to the progenitor, we find 
that these differences can reach up to 10 degrees in our experi¬ 
ments. 

In summary, stream particles are distributed following 
straight lines in angle space when computed using the correct 
potential and otherwise show a wiggly behaviour and a broken 
energy gradient. These straight lines have the same slope in fre¬ 
quency space in the static case. If the potential has evolved in 
time, then the slope in frequency space is typically shallower 
than in angle space, to an extent that depends on the growth of 
the gravitational potential. 

4.2.3. Validation of the analytic model 

We expect our analytic model to faithfully reproduce the be¬ 
haviour of the test-particle simulations except in some cases. We 
neglected two effects that can modify these slopes in the deriva¬ 
tions above. The first appears when the initial angle spread of 
the progenitor is comparable to its current extent, in which case 
the former needs to be taken into account. Since the initial extent 
typically is unknown, it is better to use very extended streams to 
derive the evolution in time of the host potential. The second ef¬ 
fect is when the time evolution of the potential is non-adiabatic. 
In this case, the higher order terms of Eq. ( p^ have to be taken 
into account. 

To confirm this, we used three simulations with different 
growth factors ag = {0.0,0.4,0.8}, but with the same initial ac¬ 
tion distribution as in the = 0.8 case by using the transforma¬ 
tion described in Eq. ( p^ . 

Figure [T^ shows the difference between the slopes in angle 
and in frequency spaces for our analytic model. In this figure, 
the green line shows the prediction of our model using Eqs. ( [25] ) 
and ( [27] ), and the actual measured slope differences from the test- 
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particle simulations are shown as blue crosses. This comparison 
shows that the agreement is excellent, as exemplified for Orbit 
07 in the left panel, and this is true for most of the orbits we 
explored. However, there are a few cases, such as those shown 
in the panels on the right, where the model fails. 

For these two cases, the initial extent of the progenitor in 
angle space cannot be neglected, as can be seen from the bottom 
panels. This can either be because the streams are relatively short 
or because the initial spread in angles is large. When the initial 
angle spreads as given by Eqs. and ( [^ are removed, the 
simulations match the model significantly better, as indicated by 
the red crosses. 

The rightmost panel still shows some disagreement, espe¬ 
cially for higher values of the growth factor ag. This behaviour 
stems from the non-adiabatic evolution of the orbit, so that 
higher terms in Eq. ( [T^ , and the last approximation in Eq. ( j^ , 
fail. An indication that these orbits evolved non-adiabatically 
is for example that Jr(0)/Jr(tf) - I ^ 0.18 for orbit 015. The 
orbits that are more circular typically suffer more strongly from 
this effect. This implies that radial orbits may be preferred to 
obtain the amount of time evolution for a stream that evolved in 
a smoothly growing potential such as ours. 

5. Observational prospects 

Now that we have determined the signature of time evolution 
of a gravitational potential through the difference in slopes in 
the angle and frequency spaces, we are interested in establishing 
whether this effect is measurable. This is particularly relevant in 
the context of the upcoming catalogues from the Gaia satellite. 

We therefore convolved our simulated data in observable 
space with errors. The largest uncertainties typically come from 
the errors in distance and proper motion. We investigated two 
cases for the errors in the parallax: 1% (cTj^ln = 0.01) and 
10% ((Tj:l7i = 0.1), while the proper motion errors are set to 
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Fig. 16: Angle and frequency space plots and residuals of the frequency space fits of all 15 streams for 1% distance (and proper motion) errors 
in the time-dependent potential. The streams are centred on the position of the centre of mass of the progenitor in the error-free case. The black 
points are error-convolved, while those in yellow are error-free and only shown in the bottom row. The residuals were computed by subtracting 
the best-fit straight line to the error-free frequency distribution (green line) from AH^. The text in the panels shows the relative difference in the 
slopes for the error-convolved and error-free cases derived in angle space Eg and in frequency space Fq (i.e. using the red vs green straight lines, 
respectively). The panels for which Et < 0.005 are marked with a green box. The difference in slope in angle space clearly is generally far below 
this value, while the frequencies are much more strongly affected by the errors. 


be cr^ = 0.5(T;r), as given by the Gaia mission error estimates. 
We assumed radial velocity errors of 1 km/s, which is the level 
obtainable from follow-up surveys of the Gaia satellite such as 


4MOST dde Jong et and WEAVE (|Dalton et al.|20T2li. 


Although this error may be seen as small, we have found no sig¬ 


nificant dependence on the radial velocity error, even if this was 
as large as 10 km/s. After applying error convolution, we con¬ 
verted the data back to Cartesian coordinates. 

In Fig. [T^ we show the structure in angle (top row) and in 
frequency (middle row) space for the 1 % distance errors for all 
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but now with 10% parallax (and proper motion) errors. The colours, lines, and insets are 


Fig. 17: Selection of the seven best orbits of Fig. 

the same as in Fig.[^ With these lOx larger distance errors, the angle slopes of the streams are still well measured for most orbits, while the 
frequency slopes generally are not. 
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Fig. 18: Selection of the orbits of Fig. [T^ with 10% parallax errors, but now with two cuts: cr^ < 2 kpc and - (Q^) | < SD^^^. The colours are 
the same as in Fig. Eli with the additionmat green points represent the particles that remain after the cuts. The grey vertical dashed lines indicate 
(Qr) and (Qr) ± applied after the distance cut. The insets in the top panels show the number of remaining particles after the cuts Ncm, while 
those in the bottom panels are the number of particles within a distance of 20 kpc from the Sun N{d < 20 kpc). 
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15 orbits evolved in the time-dependent potential. The bottom 
row panels show the residuals in frequency space after subtract¬ 
ing the fitted slope from the error-free case (green line). Here 
the black points correspond to the error-convolved case, while 
the yellow points are error free. 

We define the relative difference in slope in the frequency 
and angle Eq spaces after error convolution as 


Ea = 
Ee = 


5(AQ)-5(AQ) 

S(AD) 

S(A6)-S(A6) 
S(A0) ’ 


(36) 


with the S representing the slopes derived by fitting the distri¬ 
butions with errors. The various panels in Fig. [T^ show that this 
relative difference typically varies from 10“^ up to 6 x 10“"^ for 
Eg and between 3 x 10“"^ up to 10“^ for indicating that errors 
affect the distribution of frequencies more severely than that of 
angles. 

An estimate of tolerable errors may be derived from the fol¬ 
lowing arguments. If A^ is the slope difference between fre¬ 
quency and angle space (typically a few %), then we can tolerate 


^£2 + £2 < A5/5 « 0.005. (37) 

The panels that satisfy this requirement are highlighted with a 
green box. The very radial orbits that have small pericentres 
are mostly selected in this way. These streams have many stars 
relatively near the Sun (which is positioned at 8 kpc from the 
Galactic centre along the x axis). The more radial streams also 
have a much broader spread in which is slightly less affected 
by the errors than the more circular orbits, which have a broader 
spread in angular momentum. We conclude that the 1% distance 
errors allow us to uncover the time evolution for at least half of 
our orbits when the growth factor is = 0.8. However, if we 
were to consider a growth factor = 0.4, the misalignment be¬ 
tween angle and frequency space would typically be reduced by 
a factor 2, and the errors would considerably hinder measuring 
the time evolution for more slowly evolving potentials. 

In Fig. [T^ we show the seven best orbits from Fig. [T^ but 
now with 10% distance errors. The colours and insets are the 
same as in Fig. [T^ As the insets and green boxes show, the an¬ 
gles are still quite well measured, but the frequency slope is not. 
Only one very radial orbit (that comes close to the Sun) presents 
an acceptable level of uncertainty. There is significant structure 
in the bottom panels, which depict the residuals between the dis¬ 
tribution in frequency space before and after error convolution. 
For each experiment this is a result of superposing individual 
wraps with varying distance gradients, and this, after a distance- 
dependent error convolution, causes different frequency distri¬ 
butions. The bow-tie-like structure is thus the result of the over¬ 
lap of multiple wraps that are differently distorted in frequency 
space. The more circular orbits 015 and 09 only show one leg of 
these bow-tie-like structures because basically just slightly more 
than one wrap is present for these experiments. 

Given that the 10% errors make it challenging to fit a slope 
in frequency space, but less so in angle space, we investigated 
methods for cutting the data and retaining the highest quality 
measurements. Figure \TE\ shows the result of combining a dis¬ 
tance cut and a frequency cut for the same seven orbits as in 
Fig.[T^ The cuts are given by 

(Td <2 kpc, 

|Or-<a)l< SDa^, (38) 


with (O^) the mean frequency from the data and the stan¬ 
dard deviation of both computed after the distance cut. The 
distance error cut improves both the tangential velocities and the 
distances, and its effect is to mainly select particles at pericentre, 
as seen in the top panels of this figure. Our cut at 2 kpc essen¬ 
tially removes particles with errors larger than 10% at 20 kpc. 
The motivation for imposing a frequency cut is to remove out¬ 
liers. The vertical dashed lines in the frequency space of Fig.p^ 
show the mean and standard deviation of AQ^ with errors, and 
they encompass the error-free AO^ distribution indicated in yel¬ 
low quite well. The combination of these cuts still leaves many 
particles in the stream, as indicated by the insets in the figure, 
except for orbits 015 and 09, which basically have no or too 
few particles closer than 20 kpc, and hence have distance errors 
larger than 2 kpc. In conclusion, a cut in frequency and distance 
leads to satisfactory improvements, except when there are too 
few particles and for a few pathological cases. Such a cut might 
therefore be used to measure the slope in frequency space, but is 
not necessary in angle space. 


6. Discussion and conciusions 

We have investigated the effect of an adiabatically growing 
time-dependent potential on streams. To this end we used the 
inside-out growing spherical NFW potential from |Buist & Helmi 
( |2014| ) as a background for modelling the evolution of a set of 
15 test-particle streams. 

We performed a series of numerical experiments starting 
from different initial conditions that were mostly ran in a quite 
strongly growing model in which the enclosed mass within the 
orbits that we explored approximately doubled. These experi¬ 
ments showed that the precession rate of streams, that is, the an¬ 
gular location of the ‘apocentres’ of streams in the orbital plane, 
is significantly different when the potential has evolved in time. 
Time evolution typically leads to a misalignment or angular dif¬ 
ference of ~ 10° in comparison to the static case, and this is 
roughly independent of the progenitor size. To be able to detect 
this effect, however, streams need to be sufficiently long, mean¬ 
ing that they need to have wrapped around more than once. 

We then analysed the behaviour of streams in action-angle 
coordinates. Streams typically appear as extended linear struc¬ 
tures in angle and frequency space, also in the time-dependent 
case. We found that time evolution causes these distributions 
to differ in slope, unlike what is expected for a static potential. 
To explain these findings, we developed an analytic description 
of a stream in an adiabatically changing time-dependent poten¬ 
tial. This allowed us to predict that the difference in slopes in 
angle and in frequency space is of order ~ 0.005 - 0.025, that 
is, 0.3° - 1.5°, for potentials with a growth factor ranging from 
0.2 to 0.8. Experiments for which the time evolution signal is 
most reliably recovered were those that initially had a small size 
and evolved in an adiabatic way. This condition is most easily 
satisfied for the more radial orbits. 

Although the predicted effect is small, we explored whether 
it would be observable with the next-generation facilit ies for 
mapping the Galaxy, such as the Gaia satellite mission ( [Prusti 
|2012| ) and the follow-up surveys 4MQST ( |de Jong et al.||2012| ) 
and WEAVE ( [Dalton et al.||2012 ). When assuming an error in 
parallax of 1% (and half of that in the proper motion), we found 
that for about half of the streams the time-dependence signa¬ 
ture is expected to be measurable. When the errors in parallax 
(and consequently in proper motion) are increased to 10%, only 
1 of our 15 streams can be used to determine the time evolu¬ 
tion. Reasonable additional cuts in distance errors and in fre- 
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quency space (because this is most strongly affected by errors), 
allowed the time evolution to be determined for an additional 
four streams. The best results were obtained for streams that 
are very radial and have small pericentres (inside the solar cir¬ 
cle). This is predominantly because these streams have many 
stars close to the Sun and therefore have smaller observational 
uncertainties. There is no guarantee that streams with such char¬ 
acteristics will be present in the Gaia dataset, of course, both 
because of the requirement on their orbital properties and on the 
rather optimistic assumptions on the distance and proper motion 
errors. 

Another important issue is that the correct potential is nec¬ 
essary to derive the correct angle and frequency space distribu¬ 
tions. Fortunately, an incorrect potential distorts the distribu¬ 
tion in angle and in frequency spaces in such a manner that each 
space can be used independently to improve the estimation of 
the potential for sufficiently long streams. For short streams (e.g. 
with only one wrap), the misalignment observed might be con¬ 
fused with time evolution instead of the incorrect assumptions 
about the potential. The use of short streams might therefore in¬ 
duce a bias in the derived values of the characteristic parameters 
of the potential or its time dependence. 

By construction, we neglected the effect of self-gravity in our 
experiments. Self-gravity would create a gap in the action space 
between (the progeni tor and) the leading and trailing streams 
( Gibbons et al.||2Q14| ), because the particles released at specific 
points along the orbit (rather than continuously) are offset in en¬ 
ergy ( |Johnston||1998| ). This structure will most likely also be 
apparent in the angle and frequency distributions of the particles 
that make up the stream ( [Eyre & Binney|2QlT] ). Moreover, when 
close to the progenitor, released particles still exp erience a grav¬ 
itational pull by the progenitor ( |Choi et al.|2009| ), and epicyclic 
oscillations are seen in streams in N-body simulations ( |Kupper| 
et al.|201()||2012| ). These effects may complicate fitting a straight 


line to the distribution and determining its slope accurately, but 
after the particles become unbound, the dynamics is essentially 

the same as in our test-particle integrations (ISanders & Binneyl 

- 

Other interesting time-dependent effects may be involved 
in shaping streams and complicate interpretations. Interactions 
with dark subhaloes orbiting the Galactic halo will affect the 
structure of streams and create gaps, for exam ple ( |Yoon et ak] 
2011[ |Carlberg||2013[ |Ngan & Carlberg|[2Q14| ). Dwarf galax¬ 


ies as massive as the Large and Small Magellanic Clouds are 
able t o significantly perturb the potential ( |Vera-Ciro & Helrni] 
2Q13|), and might induce non-adiabatic changes in the stream or- 


bits ( [Gomez et al.|2015| ). It therefore seems important to try and 
understand these less secular effects on the dynamical evolution 
of streams. 

We have presented here the first steps towards understanding 
the imprints of time evolution. More explorations are necessary 
because we did not study other density profiles or deviations 
from spherical symmetry, nor did we take into account the 
effects of a five halo or a (growing) disk. It may be possible 
to extend our method to any potential for which action-angle 
coordinates can be derived or approximated, and we do expect 
the general behaviour to be similar when the potential is growing 
adiabatically. Overall, we expect the imprint of the smooth mass 
growth of the Galactic dark halo to be present in streams and to 
a have a small but non-negligible magnitude. 
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Appendix A: Computing the angies 

The angles result from the derivatives of the generating function 
W(x, J) w.r.t. the actions 


The derivative Wr,H has already been worked out for Or, and the 
other derivatives are 


6 = 


dW(xJ) 
dj ’ 


(A.1) 


•/ 


= -^\ -2=-^ 


dr _ ifpr>0. 

Pi- \2/2(?-apo) - flir) if Pr < 0, 


( |Binney & Tremaine|2008[ Eq. 3.204). The complete general- 
ing function for a spherical system is given by (modified from 
Eq. 3.220, [Binney & Tremaine|2008| ) 

W(x, J) = W^((P, J) + J) + Wr(r, J) 

pr 

= (pj(j)+ I /^) + I dr Pr(Jr, J(f), 

^^min *^^peri 

(A.2) 

where we integrate over the particles trajectory in phase space, 
which means for example that during a whole radial period, the 
particle transverses twice the branch from pericentre to apocen- 
tre, but once in reversed direction. The latitudinal momentum 
is 




Pr 


/2(?'apo) = 


(A.9) 




{-V 

J P9 V 


/ 


O O L- / COS i 

pI = L^ - ^=L^l -^ 

sin^ d \ sin^ d 

where cos i = L^jL. The radial momentum pr is 


sin^ J J ^ 

MP)= f 

/3(^max) = 0. 


-fsid) ifp^<0 


1 


1^ 1 

d§ 


COS^ I 

sin^ & 

cos i 
sin^ 


V I _ cos^ i 
sin^ 


(A.IO) 


(A.3) 


(A.4) 


p2=2(//(J)-0(r))-—, 

ri 

where L = Jg + \J^\ and 7^ > 0 is assumecj^ 

We find the radial angle as 

= (mod 2.), (A.5) 


The function W^,/^ is oscillatory in nature (w.r.t. d), while the 
combination Wr,j^ + ^cpWr^n is also oscillatory (w.r.t. r): after 
one radial period it evaluates to - 2 L/ 2 (rapo) + 2 ^ 0/1 (rapo) = 0. 
The latitudinal angle is given by 


dW(x, J) 

dh 


= + ^§Wr,H, (mod In). 


(A.11) 


dJr 


dJr OH 


Because we assume J^p > 0, we find Wr,j^ = Wr,j^ and 
The remaining derivative of the generating function is 


where we note that the angles are always defined modulo 2n 
because there is no information on how many loops a particle 
made around the Galactic centre. The derivative of the generat¬ 
ing function is given by 




cdd r, 


1 


V I _ cos^ i 


jf4(d) if p» > 0, 
127r - /4(j?) if pg < 0, 


W, 


r,H 


r dr _ 

(flir) 

if Pr > 0, 

fm= d§- 

1 

1 0 ■ 

J Pr 

l2/i(rapo)-/i(r) 

if Pr < 0, 

d&nim 1 

I COS"^ 1 


flir) 


-f 

P net 


dr 

Pr' 


.A(^max) — 


(A.12) 


/l(?-apo) = 


(A.6) 


where the conditions on pr are necessary to take the right branch 
of fi . After one full period, we find Or = 2n as expected. 

For the azimuthal angle we find 

dWix T) 

^ + ^<^Wr,H, (mod 2n), 

OJd) 


where we find 0 using the (quadrant-aware) arctangent 
0 = arctan(^/x). 


where the term Wr,j^ + ^§Wr,H vanishes after one radial period, 
while the term contains the dependence on d and increases 
by 2n after one period in d. 


Appendix B: Transformation equations 

The linearised transformation between action-angle coordinates 
and Cartesian coordinates is 


(A.7) 

T = 

■ d^W dj d^W 

1 

d^W dJ 

djdj d{\ did{\ 
dl 

dm 5p 

(A.8) 

d(\ 

5p 


^ Otherwise, every derivative of L w.r.t. Jp results in a multiplication 
with sign(/ 0 ). 


(B.l) 


For simplicity, we provide here the matrix for the 2D case (i.e. 
when the orbit is in the plane). In that case, Je = ^ and = L = 
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L^, so that we find 


T = 


1 tOcj^r ^0^,pr 

0 %,r t0^,p^ t0^,Pr 

0 0 1 0 


(B.2) 


0 tjr,r hr,Pcl, 


^Jr, Pr 


where the above terms are given by 
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The functions k and p are 



^ L 

k = £1^--^. 


(B.8) 

5<D l2 
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(B.9) 


and the W /.^ j. are found by differentiating the generating func¬ 
tion 
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dJ] Jr^^pAdJ^ r2 pA' 
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